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Abstract 



The present analytical understanding on the nature of the singularities which form 
at the endstate of gravitational collapse of massive fluid bodies ("stars") is reviewed. 
Special emphasis is devoted to the issue of physical reasonability of the models. 
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The investigation on the "final fate" of gravitational collapse of initially regular distributions 



m 
o 

Q\ [ of matter is one of the most active field of research in contemporary general relativity. It is, 
indeed, known that under fairly general hypothesis solutions of the Einstein equations with 
"physically reasonable" matter can develop into singularities The key problem that still 
remains unsolved is the nature of such singularities. The main open issue is whether the 
bJT)! singularities, which arise as the end point of collapse, can actually be observed. 

Roger Penrose @, was the first to propose the idea, known as cosmic censorship con- 
jecture: does there exist a cosmic censor who forbids the occurrence of naked singularities, 



clothing each one in an absolute event horizon? This conjecture can be formulated in a 
"strong" sense (in a "reasonable" spacetime we cannot have a naked singularity) or in a 
weak sense (even if such singularities occur they are safely hidden behind an event horizon, 
and therefore cannot communicate with far-away observers). Since Penrose's proposal there 
have been various attempts to prove the conjecture (see || and references therein). Unfor- 
tunately, no such attempts have been successful so far. As a consequence, the research in 
this field turned to more tractable objectives. In particular, one would like to understand 
what happens in simple systems, like spherically symmetric ones (interestingly enough, even 
this apparently innocuous problem is far from being completely solved, although, as we shall 
sec. a general pattern does seem to arise). 

Our aim here is to overview only models which have a clear physical interpretation. 
Therefore, we shall require satisfaction of the weak energy condition as well as existence of 
a singularity free initial data surface. Moreover, we shall take into consideration solutions of 
the Einstein field equations which are physically meaningful in terms of a (phenomenological) 
equation of state of the matter. In this respect it is worth mentioning that material schemes 
having a well defined microscopical interpretation, like the Vlasov-Einstein system, would be 
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closer to such a requisite H] . However, very little is known on the dynamics of such models 
(a numerical investigation has been recently carried out ||). 

In presence of excellent general reviews on gravitational collapse and cosmic censorship || 
[7|, §, ||, [H| we have focused our attention on a specific issue namely, the investigation of 
analytical models describing the gravitational collapse of massive stars. Therefore we are not 
going to address here many related important topics. These include Vaidya spacetimes H, 
radiation shells (see |§ and references therein), gravitational collapse of scalar fields Jl2 



13 1, critical behaviour in numerical relativity WM, stability of Cauchy horizon in Reissner- 



Nordstrom spacetimes [fl5|l, the Hoop conjecture [O, W% among others fi [TO, [T3, II 



2 Einstein equations for spherical collapse 

What is known analytically in gravitational collapse is essentially restricted to spherical 
symmetry (one exception is given by the Szekeres "quasi-spherical" spacetimes f2~0fl , but the 
results |2l| , |22|j are very similar to those holding in Tolman-Bondi models - See Section 4.2.4). 
Therefore we discuss the mathematical structure of the Einstein field equations describing 
collapse of a deformable body only in the spherically symmetric case. For perfect fluids this 



structure is well known fl2~3j; we present here a more general case which takes into account 



anisotropic materials as well |2~4" |. Also, we consider only non-dissipative processes since very 



little is known in the dissipative case p5 . 



The general, spherically symmetric, non-static line element in comoving coordinates 
t, r, 6, tp can be written in terms of three functions z/, rj, Y of r and t: 

d s 2 = -e 2v dt 2 + rj^dr 2 + Y 2 (d6 2 + sin 2 6d<p 2 ) . (1) 

Throughout this paper, we shall assume that the collapsing body is "materially spherically 
symmetric" in the sense that all the physical "observables" do not depend on angles. The 
matter density of the material (baryon number density) is then given by 

p = po(r)Y~ 2 ^ , (2) 

where po(r) is an arbitrary (positive) function. As in any theory of continuous media, to 
choose a specific material one has to specify the internal energy e. This function depends on 
the parameters characterizing the state of strain of the body (for a recent review of relativistic 



elasticity theory see [0). It can be shown that such parameters can be identified with Y and 
rj in the comoving frame (in other words, any deformation is described "gravitationally"). 
Therefore, we introduce as equation of state of the material a positive function 

e = ip(r, Y, 77) . 

Here the explicit dependence on r takes into account possible inhomogeneities. 

The energy-momentum tensor can be readily calculated and the result is a diagonal tensor 
of the form = diag(— e, S, n, n). Here the stress-strain relations (i.e. the relations giving 
the radial stress £ and the tangential stress n in terms of the constitutive function) are given 
by 

s = 2^ - i> , (3) 
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We shall always use the word stress rather than pressure since both E and IT can be in 
principle negative (tensions) without violating the energy conditions. 

Different materials correspond to different choices of the function ip. It is, however, worth 
mentioning that the material scheme most widely used in astrophysical applications is the 
perfect fluid, which can be characterized as a material whose function of state depends on 
the number density only (e = ip(p))- In this case both stresses coincide: 

dtp ~ 

n = E = p- ip :=p , 

dp 

where p is the isotropic pressure. Two particular cases are worth mentioning in this scheme. 
One is that of linear pressure-density relation (p = ke); the corresponding equation of state 
is ip(p) = p(l + Ap k ) where A is a constant. The other one is the dust model for which p = 0. 
In this case the energy is distributed proportionally to the mass (ip(p) = p). 

As soon as one allows anisotropy to occur, other interesting models appear (see J27] 



for a review on the role of anisotropy in relativistic astrophysics). Recently, a particular 



anisotropic model has been singled out |28| (for previous investigations on this kind of 
models see references in p8|). In this model, one assumes that the radial stress identically 
vanishes. The key role is played by the equation (pj) which shows that the dependence of ip 
on rj must be a multiplicative dependence from ^Jfj only. Therefore, materials with vanishing 
radial stresses can be characterized, using equation (^), via equations of state of the form 

ip(r,ri,Y) = ph(r,Y), (5) 

where h is a positive, but otherwise arbitrary, function. 

Once an equation of state has been chosen, the Einstein field equations become a closed 
system; in spherical symmetry there are three independent equations for the three variables 
v, f] and Y. It has proven, however, to be very useful to write the field equations as a system 
of four differential equations. This is done by introducing the mass function p3| , ^ defined 
as 

m (r, t) = ^-(l- Y'\ + Y 2 e~ 2u ) , (6) 
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where a dash and a dot denote derivatives with respect to r and t respectively. The mass 
function is arbitrary (positive) and allows us to write the field equations in the following 
form (four compatible equations for three variables) 

rri = AneY 2 Y' , (7) 

m = -AkYY 2 Y , (8) 

Y'r] = -2r](Y' - Yv') , (9) 

£' = — (e + £)z/ — 2(E — U)(Y'/Y) . (10) 



3 Physical reasonability and initial data 

It is easy to produce new "solutions" of the Einstein field equations in "matter". Indeed, 
just pick up a metric at will and claim that this is a "solution" referring to the calculated 
energy momentum tensor. Of course, what one has to do to remove the quotation marks 
from the above statements is to check the physical reasonability of the results. 
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First of all, the weak energy condition must be imposed on the energy momentum tensor. 
This condition requires T flu ii tJ 'u l/ > for any non spacelike w M and implies, besides positivity 
of e, non-negativeness of e + £ and e + II. Due to equations (|) and ([|), such conditions are 
equivalent to differential inequalities on the function if>, namely, |^ > and |^ < 0. 

Imposing the weak energy condition per se does not assure physical reasonability, since 
there is no guarantee that the energy momentum tensor can be deduced from a field theoretic 
description of matter. What is needed is the satisfaction of a suitable equation of state. We 
shall require, in addition, the equation of state to be locally stable [23]| (this last requirement 




could be relaxed in presence of rotationally-induced stress). The local stability condition 
requires the (local) equilibrium state of the material to be unstrained. In the "comoving 
picture" this amounts to say that the energy must have an absolute minimum at the flat- 
space values of the metric. 

Let us collect the functions describing admissible equations of state in spherical symmetry 
in a set 

Vl/ = \ip G C 2 (Rl,R + ) : if)(r,r,l) = mmip(r, Y, 77), ^ > 
[ Of] 

We shall always assume that the value of if) at the minimum (which in general can be a 
function of r) has been rescaled to unity (this can be done without loss of generality). 

A solution of the Einstein field equations describes the collapse of an initially regular 
distribution of matter only if the spacetime admits a spacelike hypersurface (t = say) 
which carries regular initial data. This means that the metric, its inverse, and the second 
fundamental form all have to be continuous at t = 0. 

On the initial hypersurface we use the scaling freedom of the r coordinate to set Y(r, 0) = 
r. We call a set of initial data complete if it is minimal, in the sense that no part of its content 
can be gauged away by a coordinate transformation. We will now prove that a complete 
set of initial data for equations (|7|)-(|10l), at fixed equation of state, is composed by a pair 
of functions. Physically, such functions describe the initial distribution of energy density 
eo = e(r, 0) and of velocity Vq = e~ u ^ r ^Y{r 1 0). It is sometimes convenient to parameterize 
these two distributions in terms of two other functions {F(r), f(r)} where F(r) = m(r, 0) 
is the initial distribution of mass and f(r) is called "energy function". The relationship 
between the two sets {F, f} and {e , Vo} is given by the following formulae 

F' = 47rr 2 e , / = V 2 - 2F/r . (11) 

The function F has to be non-negative with eo(0) = lim F(r)/r 3 finite and non- vanishing, 

while / has to be greater than —1 to preserve the signature of the metric (see equation (|13|) 
below) with lim/(r) = 0. 

r— >0 

The proof that only two arbitrary functions of r are a complete set of initial data at 
fixed equation of state is implicitly contained in many papers on spherical collapse (see 
e.g. [f23|, |3(J). It seems, however, that a complete proof has never been published in details, 
so we take this occasion to give it. We denote by a subscript, the initial value of each quantity 
appearing in the Einstein field equations. We know F = m , Vq = c~ U0 Yq. Since Y = r and 
if) is a known function of r, Y and r), its initial value ipo is also known - as well as the initial 
values of the stresses due to formulae ([|) and (^) - as a function of r and rjo (ifjo = if)(r, r, 770)). 
Evaluation of equation (|7D at t = now gives ipo = F'/(4irr 2 ), from which the value of rjo 
can be extracted algebraically. At this point, the remaining three field equations can be used 
to evaluate the remaining data, i.e., m , z/ and r) : 
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m = -47rE rV , i/ = -y^=- - 2 f ° "° , 770 = . (12) 

This completes the proof. 

In what follows, we consider only solutions which can be interpreted as models of col- 
lapsing stars, i.e. isolated objects rather than "universes". This is possible only if the 
metric matches smoothly with the Schwarzschild vacuum solution (the matching between 
two metrics is smooth if both the first and the second fundamental form are continuous on 
the matching surface). 

4 Classification and nature of singularities 

To understand the collapsing scenarios as well as the nature of the singularities, one would 
like to analyse exact solutions of the Einstein field equations. However, it goes without 
saying that the non-linearity of such equations makes them essentially untractable, even 
in spherical symmetry, without further simplifying hypotheses (like e.g. vanishing shear or 
acceleration, see [[H, 55]). In view of such difficulties, one would like to extract information 



from the Einstein equations without solving them completely. 

The mathematical structure of spherical collapse discussed in the previous Section shows 
that there is a one-to-one correspondence between solutions and choices of triplets (s, say) of 
functions s = {F(r), f(r),ip(r,r],Y)} (we are, of course, identifying solutions modulo gauge 
transformations). We denote the set of solutions parameterized in this way by S. Since 
in this parameterization we have already taken into account regularity as well as physical 
admissibility, the whole physical content of the cosmic censorship problem can be translated 
in the mathematical terms of predicting the endstate of any choice of s G S. 

4.1 Non-singular solutions. 

The existence of non-singular, non-static solutions of the Einstein field equations is not 
forbidden by the singularity theorems (we refer the reader to for a recent review on sin- 
gularity theorems and further references) and a gravitational collapse can lead to a bouncing 
back, at a finite area radius, without singularity formation. This phenomenon can be recur- 
sive, producing an eternally oscillating, globally regular solution |34], [36[. Speaking very 
roughly, one can prepare regular initial data in such a way that the region of possible trapped 
surface formation is disconnected from the data. In this way the remaining hypothesis of 
the singularity theorems can be satisfied without singularity formation (see also p4| ). 



Much more exotic than the oscillating solutions, are the globally regular solutions that 
describe non-singular blackholes f33fl . These are matter objects "sitting" inside their Schwarz- 



schild radii. Their finite extension replaces the central singularity with a matter-filled, non 
singular region. In this case trapped surfaces are obviously present and therefore the strong 
energy condition must be violated (a simple argument shows that also the dominant energy 
condition is necessarily violated). In any case, such strange objects are not ruled out by 
general relativity as far as only the weak energy condition and an equation of state are 
required |J7|, |38|, [J9|. However, as far as we are aware, it is not known if a fully dynamic 
solution exists which could eventually lead to this exotic end state. 
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4.2 Singular solutions 



We now move to the case in which a singularity is formed in the future of a regular initial data 
set. The singularities of spherically symmetric matter filled spacetimes can be recognized 
from divergence of the energy density and curvature scalars, like e.g. the Kretschmann scalar 
Jl^po 'H Essentially, these singularities can be of two kinds. We shall call shell crossing 
singularities those at which Y' vanishes (Y ^ 0), and shell focusing singularities those at 
which Y vanishes. To the two kind of singularities correspond two curves t sc {r) and £/(r) in 
the r-t plane, defined by Y'(r,t sc (r)) = and by Y(r,tf(r)) = respectively. Physically, the 
shell crossing curve gives the time at which two neighbouring shells of matter intersect each 
other, while the shell focusing curve identifies at which time the shell labeled r "crushes to 
zero size" . Of these two kinds of singularities, the one of physical interest is obviously that 
occurring first in the sense that, at fixed r, one has tf > t sc or vice versa. It would be very 
interesting, therefore, to carry out a study of the field equations in order to obtain conditions 
for shell crossing avoidance in spherical spacetimes. This study has been up to now carried 
out only for dust spacetimes |H| f|T], p£2f . 

From the point of view of censorship, the nature of a singularity in an asymptotically flat, 
initially regular spacetime can be one of the following. First of all, a singularity can be space- 
like, like e.g. the Schwarzschild singularity or the singularity occurring at the endstate of the 
collapse of a Oppenheimer-Snyder dust cloud (see Section 4.2.2 below). These singularities 
lie in the future of all possible observers and therefore are strongly censored (i.e. allowed 
by strong cosmic censorship). If a singularity is not strongly censored, then it is naked, i.e. 
visible to some observer. However, two different cases can occur, namely the singularity can 
be locally or globally naked. A singularity is locally naked if light signals can emerge from it 
but fall back without reaching any asymptotic observer. A singularity of this kind will be 
visible only to observers who have crossed the horizon; therefore, the weak cosmic censorship 
holds for such endstates. An example of a spacetime containing a locally naked singularity 
is provided by the Kerr spacetime with mass greater than the angular momentum per unit 
mass (or by the Reissner-Nordstrom spacetime with mass greater than charge). Finally, a 
singularity is globally naked if light-rays emerging from it can reach an asymptotic observer. 



4.2.1 Shell crossing singularities. 



The first explicit example showing formation of a naked singularity was found as a shell 
crossing singularity in a spherical dust cloud P5 . It can be shown that these singularities 
are timelike and are always locally naked. 

Some definitions have been proposed to put singularities "in the order of increasing 
Essentially, what is done is to check the behaviour of the invariants of 
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seriousness 

the Riemann tensor in the approach to the singularity. According to such criteria, the shell 
crossing singularities turn out to be "weak" at least as compared with the shell focusing 
singularities. This "weakness' is considered by some authors as an hint of a possible extension 
of the spacetime ^5|. However, in spite of their "weakness", there is at the moment no 
available general proof of extendibility of spacetimes through a shell cross (although some 
encouraging results exist, see pp. The unique exception is a paper by Papapetrou and 
Hamoui [47]. In this paper the authors claim to have explicitly found the extension in the 
case of "degenerate" shell crossing singularities, i.e. when the curve t sc (r) degenerates to an 
"instant of time" t sc = T = const. In this case, it is easy to check that the crossing happens 
at a "point" r = const, rather than at a "3-space" and this is the key to their treatment. 
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However, some results of this paper are unclear from the physical point of view. 

What is actually available is only a continuous extension of shell crossing singularities 
exists in the dust case [JI]. We shall show this in a slightly more general case. Integrating 



equation (H) formally with respect to time we can write 



y'2 

v- 1 = , (13) 



where / is the energy function and f2(r, t) = exp (— /q ^-dtj . Changing variable from r to 
Y = Y(r, t) one gets the metric 



ds 2 



n 2 



i + f 



Y 2 - (1 + f)n- 2 e 2u ) dt 2 - TYdYdt + dY 2 + Y 2 (d9 2 + sin 2 Odip 



If Q is finite and non- vanishing at t — t sc (r), the above metric is continuous with continuous 
inverse at such surface. This happens if v' vanishes (dust case) or, more generally, if v' goes 
to zero at least as Y' at the shell cross (this happens, for instance, in the case of vanishing 
radial stresses). 



In a recent paper [48 , Szekeres and Lun have shown that there exist a system of coordi- 



nates in which the metric is of class C 1 as t — > t . However, again, this result per se does 



not show extendibility of the spacetime (see also |49||). 



4.2.2 Non-central shell focusing singularities 

It is important to distinguish the central shell focusing singularity, i.e. that occurring at 
r = 0, from the other focusing singularities, since in many cases it is easy to prove that non 
central singularities are censored. 

A necessary condition for the visibility of a "point" r is that the condition 1 — 2 y^t) > 0> 
implying absence of trapped surfaces (see next Sub-section), is satisfied. Since the appar- 
ent horizon is the boundary of the region containing trapped surfaces, the above condition 
implies that the time of formation of the apparent horizon t a h(r), defined as 2m(r,t a f l (r)) = 
Y{r,t ah {r)), must not be before tf(r), i.e. t a h{r) > tf(r). Now suppose m(r,tf(r)) to be 
different from zero (as we have seen, this is not the case at the central singularity where m 
has to vanish as r 3 ). Then 1 — 2m(r,t)/Y(r,t) goes to minus infinity as t tends to tf(r) so 
that the singularity is covered. This shows that all the naked shell focusing singularities are 
necessarily massless, in the sense that m has to vanish there [50|. It follows that any non 



central singularity will certainly be censored if the mass is an increasing function with respect 
to t. Now equation (^j) gives rh = —4ttHY 2 Y and, since Y < during collapse, we conclude 
that non central singularities are always covered if the radial stress S is non-negative |5T| 
(in particular, all non-central singularities occurring in dust as well as in models with only 
tangential stress are covered, since the mass does not depend on time in this case). In the 
presence of radial tensions, the question is still open. It is known that a perfect fluid with 
p = ke exhibits naked non-central singularities if k < —1/3 |5Tf . 

4.2.3 Central singularities: the root equation 

The first explicit examples of the formation of a naked shell focusing singularity were pro- 
vided by Eardley [f)2|| , Eardley and Smarr |53|] and by Christodoulou f54|. Since then the 



techniques to study the nature of the singularities in spherically symmetric spacetimes have 
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been developed by many authors (see references in H3] ) and finally settled up by Dwivedi 
and Joshi |55| . 

The key idea is the following: if the singularity is visible, at least locally, there must 
exist light signals coming out from it. Therefore, by investigating the behaviour of radial 
null geodesies near the singularity, one can try to find out if outgoing null curves meet the 
singularity in the past. On such radial null geodesies, the derivative of Y(r, t) reads 

*L = Y' + r l - x / 2 Ye- v . 
dr 



Using (H) and (|13D in the above equation we obtain 



d ^- = Y> 
dr 



1 




tt 2 /2m 



/2m 
1 + 



(14) 



where Y' has to be understood as a known function of Y and r. If the singularity is naked, 
equation ( |H] ) must have at least one solution with definite, outgoing tangent at r = 0, i.e. 
a solution of the kind Y = X r a where a > 1 and X is a positive constant. Clearly, this 
behaviour is possible only if the necessary condition 1 — 2m /Y > is satisfied. Indeed Y' is 
equal to one, and therefore positive, on the initial data surface. If no shell crossing occurs, 
it remains positive, so that the right hand side of equation ([14]) cannot remain positive if 
2m I Y — 1 changes sign. 

Once the necessary condition is satisfied, one has to check if both X and a exist such 
that the solution of equation (|n|) is of the specified form near the singularity. On using 
L'Hospital rule we have 

X„ = lim (-irr§) • (15) 

Y=X r a 



+0 \ ar 



dr 



Using again (|14"D, this equation becomes an algebraic relation for X at fixed a. If a positive 
definite Xq exists the singularity is naked. 



4.2.4 The dust case 



The general exact solution of the Einstein field equations is known in the most simple case 
of vanishing pressure (dust) [56], |57|. In this case, from equation fllT)]), one gets v = 
(more precisely, e u is an arbitrary function of t only which can be rescaled to unity without 
loss of generality). Then Q = and it follows from equation (|13|) that rj" 1 = Y' 2 /(l + /). 
The mass is constant in time (m = F(r)) due to equation (Q) with £ = 0. Therefore (|6|) can 
be written as a Kepler-like equation (Y 2 — f + 2F/Y), which is integrable in parametric 
form for / ^ and in closed form for / = 0. Finally, the density can be read off from (|?]) as 
e = F'/(A7rY 2 Y'). 

A great effort has been paid to understand the nature of the central singularity in this 
solution [|4], [42], |53], |5^, ^8], |59], |6(], [H], ^2] and we now know the complete spectrum of 
possible endstates of the dust evolution in dependence of the initial data. We recall here 
what happens in the case of marginally bound solutions (/ = 0) p2] since it is sufficiently 



general to illustrate a tendency and simple enough to be recalled in a few lines. 

For marginally bound dust, the solutions s = {F, 1, 1} can be uniquely characterized by 
the expansion of the function F(r) at r = or, and that is the same, by the expansion 
of eo = F'/47rr 2 . Using this expansion in the root equation, it is not difficult to check the 
following results: 
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r (a) r (b) r (c) 

Figure 1: Penrose diagrams for collapsing dust clouds. The dotted line represents the center of 
the star, and the bold jagged line represents the singularity. The continuous line connecting i~ 
to the singularity represents the boundary of the collapsing cloud. Three different situations can 
occur, depending on the choice of initial data: (a) "generalized" Oppenheimer-Snyder collapse: the 
singularity is covered, (b) locally naked singularity (c) globally naked singularity. 

• If the first non-vanishing term corresponds to n = 1 or n = 2 equation flT5|) always has 
a real positive root: the singularity is naked; 

• If the first non-vanishing term is n = 3 the root equation reads 

2x 4 + x 3 + £x - f = , (16) 

where x 2 = X and £ = Fs/Fq 5 ^ 2 . From the theory of quartic one can show that this 
equation admits a real positive root if £ < £ c = —(26 + 15\/3)/2. Therefore, £ c is a 
critical parameter: at £ = £ c a "phase transition" occurs and the endstate of collapse 
turns from a naked singularity to a blackhole. 



• If n > 3 the limit in equation (13)) diverges: the singularity is covered. In particular, this 
case contains the solution first discovered by Oppenheimer and Snyder ]5£| describing 
a homogeneous dust cloud. 

The naked singularities mentioned above are locally naked. It can, however, be shown 
that if locally naked singularities occur in dust spacetimes, then spacetimes containing glob- 
ally naked singularities can be build up from these by matching procedures. The Penrose 
diagrams corresponding to the three different cases are shown in Figure 1. 

The above results can be extended to the general case of collapsing dust clouds, so that 
the final fate of the dust solutions s = {F, f, 1} is completely known. The final fate depends 
on a parameter which is a combination of coefficients of the expansions of F and / near r = 0, 
and a structure similar to that of marginally bound collapse arises (see [^] for details). 



4.2.5 Vanishing radial stresses 

Recently, the general solution for spherically symmetric dust has been extended to the case 
in which only the radial stress vanishes |23|, ^| . The solution can be reduced to quadratures 
using a system of coordinates first introduced by Ori |63| for charged dust. One of the new 
coordinate is the mass m which is constant in time (due to (jS|) with S = 0), the other 
coordinate is the "area radius" Y. In such coordinates the metric reads 
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ds 2 = -T 2 ( 1 - — ) dm 2 + 2Jl + f—dYdm - ^-dY 2 + Y 2 (d6 2 + sin 2 Od<p 2 ) , (17) 
Y J hu u 



where 



2m 1 + / „ . , r h du 
u = -\ -1 1 ' ■ 1 ' 



Y h 2 ' yv ' J u 2 y/TTJdm 

and the function g(m) is arbitrary. 

The problem of understanding the nature of the singularities for such solutions is essen- 
tially still open. It is, indeed, possible to write the root equation in explicit form, but this 
equation contains a sort of "non-locality" due to the integral entering in the definition of T. 
As a result only a few special cases have been analysed so far [^, |64], |65 . 



Among the solutions with tangential stresses a particularly interesting one is the Einstein 
cluster p!| |68|. This is a spherically symmetric cluster of rotating particles. The motion 
of the particles is sustained by the angular momentum L whose average effect is to introduce a 
non vanishing tangential stress in the energy-momentum tensor. The corresponding equation 



of state has the form h(m, Y) — Jl + L 2 (m)/Y 2 . Therefore, a solution is uniquely identified 
by the choice of three arbitrary functions of m only, namely F, f and L 2 (for L = one 
recovers dust). It turns out that the final state "at fixed dust background" (i.e. for fixed 
F, /) depends on the expansion of L 2 near m = (L ~ (5m y , say) |69|, [7(|. Considering, for 



simplicity, the marginally bound case, one finds that for 4/3 < y < 2 either the singularity 
does not form (the system bounces back) or a blackhole is formed. The threshold of naked 
singularity formation lies at y — 2, where a 2-parameter structure very similar to that of 
dust occurs. At y — 7/3 a sort of transition takes place and the evolution of the model is 
such that only the critical branch is changed, un-covering a part of the blackhole region in 
the corresponding dust spacetime; the non-critical branch is the same as in dust spacetimes. 
Finally for y > 7/3 the evolution always leads to the same end state of the corresponding 
dust solution. 



4.2.6 Self-similar collapse 



A spherically symmetric spacetime is self-similar if it admits an homothetic vector £, i.e. a 
vector satisfying C^g^ v = 2g fiu . In the comoving frame the dimensionless variables u, A and 
Y/r depend only on the "similarity variable" z := r/t, and the Einstein field equations be- 
come ordinary differential equations (we refer the reader to [71| for a complete treatment of 
self-similar solutions). Being governed by ordinary differential equations, self-similar spheri- 
cal collapse can be analyzed with the powerful techniques of dynamical systems theory [[72[ . 
The analysis of the singularities forming in self-similar spacetimes has been done by many 

I, rg, |75i m w 



authors for different equations of state, like dust |52 



fluids |75|, |78], |79|, radiation (Vaidya) shells |30|, ^T], Q and in general cases |83 



barotropic perfect 
541 . 



The picture arising resembles the dust case in the sense that both naked singularities and 
blackholes can form depending on the values of the parameters characterizing the solution. 



A thoroughly review of these and other features of self-similar solutions can be found in [[73] ; 
here we limit ourself to stress that naked singularities exist in self-similar solutions with 
pressure, thereby showing that the phenomenon of naked singularities formation cannot be 
considered as an artifact of dust (i.e. vanishing stresses) solutions. 
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4.2.7 General stresses 



The problem of predicting the final fate of an initially regular distribution of matter sup- 
ported by an arbitrary distribution of stress-energy (including e.g. the case of isotropic 
perfect fluids but also anisotropic crystalline structures which are thought to form at ex- 
tremely high densities), is still open even in the spherically symmetric case. First of all, one 
has to take into account the fact that there is a high degree of uncertainty in the properties 
of the equation of state at very high densities. Recently, Christodoulou |35| initiated the 
analysis of a simple model composed by a dust ("soft") phase for energy density below a 
certain value e and a stiff ("hard") phase for e > e (in the hard phase the pressure is given 
by V = e-e). 

Although the details of the collapse in presence of general matter fields are still largely 
unknown, it is very unlikely that the "embarrassing" examples of naked singularity formation 
like those occurring in dust can eventually be eliminated with the "addition" of stresses. 
Indeed, it is reasonable to think that a sector of naked singularities exists in the choice of 
initial data for any fixed equation of state The main issues that have, therefore, to be 
addressed are the genericity and the stability of naked singularity formation. 

Both the above italicized terms have a somewhat intuitive meaning that is, however, 
difficult to express in mathematical terms. Regarding genericity, one can mean that the set 
of initial data leading to naked singularities is not of measure zero. For instance, it has been 
shown |42| that among the dust solutions s = {F, f, 1} naked singularities are generic in the 
sense that at a fixed density profile F, one can always choose energy functions / leading to 
blackholes or naked singularities. A generalization of such a result would be that the naked 
singularities are generic - in this specific sense - at fixed, but arbitrary, equation of state ip 
(there is some convincing evidence for this, see [CT). 

The issue of stability is even more delicate than that of genericity. Indeed, any exact 
solution of a physical theory must survive to small but arbitrary perturbations in order to 
serve as a candidate for describing nature. In the case of exact solutions of the Einstein 
field equations, the notion of stability is very delicate due to to the gauge invariance of the 
theory with respect to spacetime diffeomorphisms. Recently, some evidence of stability of 
dust naked singularities against perturbations has been obtained [86] . 



5 Discussion: a picnic on the side of the road 

It is well known that, if the mass of a collapsing object does not fall below the neutron 
star limit (~ 3M Q ), no physical process is able to produce enough pressure to balance the 
gravitational pull so that continued gravitational collapse must occur. It is widely believed 
that the final state of this process is a blackhole. However, what general relativity actually 
predicts, in the cases which have been analysed so far, is that either a blackhole or a naked 
singularity is formed, depending on the initial distribution of density and velocity and on 
the constitutive nature of the collapsing matter. One may raise the objection that most of 
the known analytical results could be an artifact of spherical symmetry. However, from the 
numerical point of view some evidence that this is not the case is coming up. Therefore, as 
singularity theorems showed that the singularities occurring in collapse are generic and not 
any artifact of symmetry, a similar situation may hold for the nature of singularities as well. 

One may at this point ask if and when, a cosmic censorship theorem holds in nature. 
An answer to such a question could be in the negative. However, it remains to understand 
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the physics underlying the end state of gravitational collapse with respect to the choice of 
initial data at fixed matter model. In a famous novel [|87| (that inspired the film Stalker 
by A. Tarkovsky) a short visit of extraterrestrial life on the earth occurs. The gap between 
the two civilizations is so high that human beings are, with respect to the "garbage" left 
by the visitors, like ants exploring what remains on the side of the road after a human 
picnic. Something they find is useful, something useless, something even dangerous, but 
anyway everything is obscure and difficult to understand, looking like the weak shadow of a 
wonderful abyss of knowledge. Our present cosmic censorship understanding resembles this 
situation. Indeed, we are getting a variety of mathematical hints with somewhat obscure 
physical meanings. For instance: the condition on £ recalled in Section 4.2.4, the constraints 
arising in the gravitational collapse of Einstein clusters, the dimensionless numbers arising in 
Choptuik's numerical results [[14]], the condition on the radiation flux which arises in Vaidya 
collapse ||. Such "numbers" should presumably be the remnant weak shadow of a general 
theorem when its hypotheses are enormously restricted by the choice of the equation of state 
and of the adopted symmetries. 

To get rid of this puzzle appears to be one of the most exciting objectives of future 
research in classical relativity theory. 
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